E 



33 



Search for Hadronic Decays of W and Z Bosons in Photon Events in pp Collisions at 

v/i=1.96 GeV 

T. Aaltonen, 23 J. Adelman, 13 T. Akimoto, 54 M.G. Albrow, 17 B. Alvarez Gonzalez, 11 S. Amerio, 42 D. Amidei, 34 
A. Anastassov, 51 A. Annovi, 19 J. Antos, 14 M. Aoki, 24 G. Apollinari, 17 A. Apresyan, 47 T. Arisawa, 56 A. Artikov, 15 
W. Ashmanskas, 17 A. Attal, 3 A. Aurisano, 52 F. Azfar, 41 P. Azzi-Bacchetta 42 P. Azzurri, 45 N. Bacchetta, 42 
W. Badgett, 17 A. Barbaro-Galtieri, 28 V.E. Barnes, 47 B.A. Barnett, 25 S. Baroiant, 7 V. Bartsch, 30 G. Bauer, 32 

P.-H. Beauchemin, 33 F. Bedeschi, 45 P. Bednar, 14 S. Behari, 25 G. Bellettini, 45 J. Bellinger, 58 A. Belloni, 22 
D. Benjamin, 16 A. Beretvas, 17 J. Beringer, 28 T. Berry, 29 A. Bhatti, 49 M. Binkley, 17 D. Bisello, 42 I. Bizjak, 30 
R.E. Blair, 2 C. Blocker, 6 B. Blumenfeld, 25 A. Bocci, 16 A. Bodek, 48 V. Boisvert, 48 G. Bolla, 47 A. Bolshov, 32 

D. Bortoletto, 47 J. Boudreau, 46 A. Boveia, 10 B. Brau, 10 A. Bridgeman, 24 L. Brigliadori, 5 C. Bromberg, 35 
3. Brubaker, 13 J. Budagov, 15 H.S. Budd, 48 S. Budd, 24 K. Burkett, 17 G. Busetto, 42 P. Bussey, 21 A. Buzatu, 3 ' 
K. L. Byrum, 2 S. Cabrera r , 16 M. Campanelli, 35 M. Campbell, 34 F. Canelli, 17 A. Canepa, 44 D. Carlsmith, 58 
R. Carosi, 45 S. CarrihV, 18 S. Carron, 33 B. Casal, 11 M. Casarsa, 17 A. Castro, 5 P. Catastini, 45 D. Cauz, 53 

M. Cavalli-Sforza, 3 A. Cerri, 28 L. Cerrito p , 30 S.H. Chang, 27 Y.C. Chen, 1 M. Chertok, 7 G. Chiarelli, 45 
G. Chlachidzc, 17 F. Chlcbana, 17 K. Cho, 27 D. Chokheli, 15 J.P. Chou, 22 G. Choudalakis, 32 S.H. Chuang, 51 
K. Chung, 12 W.H. Chung, 58 Y.S. Chung, 48 C.I. Ciobanu, 24 M.A. Ciocci, 45 A. Clark, 20 D. Clark, 6 G. Compostella, 42 
M.E. Convery, 17 J. Conway, 7 B. Cooper, 30 K. Copic, 34 M. Cordelli, 19 G. Cortiana, 42 F. Crescioli, 45 
C. Cuenca Almenar r , 7 J. Cuevas , 11 R. Culbertson, 17 J.C. Cully, 34 D. Dagenhart, 17 M. Datta, 17 T. Davies, 21 
P. de Barbara, 48 S. De Cecco, 50 A. Deisher, 28 G. De Lentdecker d , 48 G. De Lorenzo, 3 M. Dell'Orso, 45 L. Demortier, 49 
J. Deng, 16 M. Deninno, 5 D. De Pedis, 50 P.F. Derwent, 17 CP. Di Giovanni, 43 C. Dionisi, 50 B. Di Ruzza, 53 
J.R. Dittmann, 4 M. D'Onofrio, 3 S. Donati, 45 P. Dong, 8 J. Donini, 42 T. Dorigo, 42 S. Dube, 51 J. Efron, 38 
R. Erbacher, 7 D. Errede, 24 S. Errede, 24 R. Eusebi, 17 H.C. Fang, 28 S. Farrington, 29 W.T. Fedorko, 13 R.G. Feild, 59 
M. Feindt, 26 J.P. Fernandez, 31 C. Ferrazza, 45 R. Field, 18 G. Flanagan, 47 R. Forrest, 7 S. Forrester, 7 M. Franklin, 22 
J.C. Freeman, 28 I. Furic, 18 M. Gallinaro, 49 J. Galyardt, 12 F. Garberson, 10 J.E. Garcia, 45 A.F. Garfmkel, 47 
K. Genser, 17 H. Gerberich, 24 D. Gerdes, 34 S. Giagu, 50 V. Giakoumopolou a , 45 P. Giannetti, 45 K. Gibson, 46 
J.L. Gimmell, 48 CM. Ginsburg, 17 N. Giokaris , 15 M. Giordani, 53 P. Giromini, 19 M. Giunta, 45 V. Glagolev, 15 

D. Glcnzinski, 17 M. Gold, 36 N. Goldschmidt, 18 A. Golossanov, 17 G. Gomez, 11 G. Gomez- Ceballos, 32 
M. Goncharov, 52 O. Gonzalez, 31 I. Gorelov, 36 A.T. Goshaw, 16 K. Goulianos, 49 A. Gresele, 42 S. Grinstein, 22 
C. Grosso-Pilcher, 13 R.C Group, 17 U. Grundler, 24 J. Guimaraes da Costa, 22 Z. Gunay-Unalan, 35 C. Haber, 28 
K. Hahn, 32 S.R. Hahn, 17 E. Halkiadakis, 51 A. Hamilton, 20 B.-Y. Han, 48 JY. Han, 48 R. Handler, 58 F. Happacher, 19 
K. Hara, 54 D. Hare, 51 M. Hare, 55 S. Harper, 41 R.F. Harr, 57 R.M. Harris, 17 M. Hartz, 46 K. Hatakeyama, 49 
J. Hauser, 8 C. Hays, 41 M. Heck, 26 A. Heijboer 44 B. Heinemann, 28 J. Heinrich, 44 C. Henderson, 32 M. Herndon, 58 
J. Heuser, 26 S. Hewamanage, 4 D. Hidas, 16 CS. Hill c , 10 D. Hirschbuchl, 26 A. Hocker, 17 S. Hou, 1 M. Houlden, 29 
S.-C. Hsu, 9 B.T. Huffman, 41 R.E. Hughes, 38 U. Husemann, 59 J. Huston, 35 J. Incandela, 10 G. Introzzi, 45 M. Iori, 50 

A. Ivanov, 7 B. Iyutin, 32 E. James, 17 B. Jayatilaka, 16 D. Jeans, 50 E.J. Jeon, 27 S. Jindariani, 18 W. Johnson, 7 
M. Jones, 47 K.K. Joo, 27 S.Y. Jun, 12 J.E. Jung, 27 T.R. Junk, 24 T. Kamon, 52 D. Kar, 18 P.E. Karchin, 57 Y. Kato, 40 
R. Kephart, 17 U. Kerzel, 26 V. Khotilovich, 52 B. Kilminster, 38 D.H. Kim, 27 H.S. Kim, 27 J.E. Kim, 27 M.J. Kim, 17 

S.B. Kim, 27 S.H. Kim, 54 Y.K. Kim, 13 N. Kimura, 54 L. Kirsch, 6 S. Klimenko, 18 M. Klute, 32 B. Knuteson, 32 
B.R. Ko, 16 S.A. Koay, 10 K. Kondo, 56 D.J. Kong, 27 J. Konigsberg, 18 A. Korytov, 18 A.V. Kotwal, 16 J. Kraus, 24 
M. Kreps, 26 J. Kroll, 44 N. Krumnack, 4 M. Kruse, 16 V. Krutelyov, 10 T. Kubo, 54 S. E. Kuhlmann, 2 T. Kuhr, 26 
N.P. Kulkarni, 57 Y. Kusakabe, 56 S. Kwang, 13 A.T. Laasanen, 47 S. Lai, 33 S. Lami, 45 S. Lammel, 17 M. Lancaster, 30 
R.L. Lander, 7 K. Lannon, 38 A. Lath, 51 G. Latino, 45 I. Lazzizzera, 42 T. LeCompte, 2 J. Lee, 48 J. Lee, 27 Y.J. Lee, 27 
S.W. Lee 9 , 52 R. Lefevre, 20 N. Leonardo, 32 S. Leone, 45 S. Levy, 13 J.D. Lewis, 17 C. Lin, 59 CS. Lin, 28 J. Linacre, 41 
M. Lindgren, 17 E. Lipeles, 9 A. Lister, 7 D.O. Litvintsev, 17 T. Liu, 17 N.S. Lockyer, 44 A. Loginov, 59 M. Loreti, 42 
L. Lovas, 14 R.-S. Lu, 1 D. Lucchesi, 42 J. Lueck, 26 C. Luci, 50 P. Lujan, 28 P. Lukens, 17 G. Lungu, 18 L. Lyons, 41 
J. Lys, 28 R. Lysak, 14 E. Lytken, 47 P. Mack, 26 D. MacQueen, 33 R. Madrak, 17 K. Maeshima, 17 K. Makhoul, 32 
T. Maki, 23 P. Maksimovic, 25 S. Malde, 41 S. Malik, 30 G. Manca, 29 A. Manousakis , 15 F. Margaroli, 47 C. Marino, 26 
CP. Marino, 24 A. Martin, 59 M. Martin, 25 V. Martin- 7 ', 21 M. Martinez, 3 R. Martmez-Ballarin, 31 T. Maruyama, 54 
P. Mastrandrea, 50 T. Masubuchi, 54 M.E. Mattson, 57 P. Mazzanti, 5 K.S. McFarland, 48 P. Mclntyre, 52 R. McNulty\ 29 
A. Mehta, 29 P. Mehtala, 23 S. Menzemer fc , n A. Menzione, 45 P. Merkel, 47 C. Mesropian, 49 A. Messina, 35 T. Miao, 17 
N. Miladinovic, 6 J. Miles, 32 R. Miller, 35 C. Mills, 22 M. Milnik, 26 A. Mitra, 1 G. Mitselmakher, 18 H. Miyake, 54 
S. Moed, 22 N. Moggi, 5 CS. Moon, 27 R. Moore, 17 M.J. Morello, 45 P. Movilla Fernandez, 28 J. Mulmenstadt, 28 



2 



A. Mukherjee, 17 Th. Muller, 26 R. Mumford, 25 P. Murat, 17 M. Mussini, 5 J. Nachtman, 17 Y. Nagai, 54 A. Nagano, 54 
J. Naganoma, 56 K. Nakamura, 54 I. Nakano, 39 A. Napier, 55 V. Necula, 16 C. Neu, 44 M.S. Neubauer, 24 J. Nielsen^, 28 

L. Nodulman, 2 M. Norman, 9 O. Norniclla, 24 E. Nurse, 30 S.H. Oh, 16 Y.D. Oh, 27 I. Oksuzian, 18 T. Okusawa, 40 
R. Oldeman, 29 R. Orava, 23 K. Osterberg, 23 S. Pagan Griso, 42 C. Pagliarone, 45 E. Palcncia, 17 V. Papadimitriou, 17 
A. Papaikonomou, 26 A. A. Paramonov, 13 B. Parks, 38 S. Pashapour, 33 J. Patrick, 17 G. Pauletta, 53 M. Paulini, 12 
C. Paus, 32 D.E. Pcllctt, 7 A. Pcnzo, 53 T.J. Phillips, 16 G. Piacentino, 45 J. Picdra, 43 L. Pinera, 18 K. Pitts, 24 
C. Plagcr, 8 L. Pondrom, 58 X. Portell, 3 O. Poukhov, 15 N. Pounder, 41 F. Prakoshyn, 15 A. Pronko, 17 J. Proudfoot, 2 
F. PtohosY 7 G. Punzi, 45 J. Purslcy, 58 J. Rademacker c , 41 A. Rahaman, 46 V. Ramakrishnan, 58 N. Ranjan, 47 
I. Redondo, 31 B. Reisert, 17 V. Rekovic, 36 P. Ronton, 41 M. Rescigno, 50 S. Richter, 26 F. Rimondi, 5 L. Ristori, 45 
A. Robson, 21 T. Rodrigo, 11 E. Rogers, 24 S. Rolli, 55 R. Roser, 17 M. Rossi, 53 R. Rossin, 10 P. Roy, 33 A. Ruiz, 11 
J. Russ, 12 V. Rusu, 17 H. Saarikko, 23 A. Safonov, 52 W.K. Sakumoto, 48 G. Salamanna, 50 O. Salto, 3 L. Santi, 53 

S. Sarkar, 50 L. Sartori, 45 K. Sato, 17 A. Savoy-Navarro, 43 T. Scheidlc, 26 P. Schlabach, 17 E.E. Schmidt, 17 
M.A. Schmidt, 13 M.P. Schmidt*, 59 M. Schmitt, 37 T. Schwarz, 7 L. Scodellaro, 11 A.L. Scott, 10 A. Scribano, 45 
F. Scuri, 45 A. Sedov, 47 S. Seidel, 36 Y. Seiya, 40 A. Semenov, 15 L. Sexton-Kcnncdy, 17 A. Sfyrla, 20 S.Z. Shalhout, 57 

M.D. Shapiro, 28 T. Shears, 29 P.F. Shepard, 46 D. Sherman, 22 M. Shimojima™, 54 M. Shochct, 13 Y. Shon, 58 
I. Shreyber, 20 A. Sidoti, 45 P. Sincrvo, 33 A. Sisakyan, 15 A.J. Slaughter, 17 J. Slaunwhite, 38 K. Sliwa, 55 J.R. Smith, 7 
F.D. Snider, 17 R. Snihur, 33 M. Soderberg, 34 A. Soha, 7 S. Somalwar, 51 V. Sorin, 35 J. Spalding, 17 F. Spinella, 45 
T. Spreitzer, 33 P. Squillacioti, 45 M. Stanitzki, 59 R. St. Denis, 21 B. Stelzer, 8 O. Stelzer-Chilton, 41 D. Stcntz, 37 
J. Strologas, 36 D. Stuart, 10 J.S. Suh, 27 A. Sukhanov, 18 H. Sun, 55 I. Suslov, 15 T. Suzuki, 54 A. Taffard 6 , 24 
R. Takashima, 39 Y. Takeuchi, 54 R. Tanaka, 39 M. Tecchio, 34 P.K. Teng, 1 K. Terashi, 49 J. Thom^, 17 
A.S. Thompson, 21 G.A. Thompson, 24 E. Thomson, 44 P. Tipton, 59 V. Tiwari, 12 S. Tkaczyk, 17 D. Toback, 52 
S. Tokar, 14 K. Tollefson, 35 T. Tomura, 54 D. Tonelli, 17 S. Torre, 19 D. Torretta, 17 S. Tourneur, 43 W. Trischuk, 33 
Y. Tu, 44 N. Turini, 45 F. Ukegawa, 54 S. Uozumi, 54 S. Vallecorsa, 20 N. van Remortcl, 23 A. Varganov, 34 
E. Vataga, 36 F. Vazquez', 18 G. Velev, 17 C. Vellidis", 45 V. Veszpremi, 47 M. Vidal, 31 R. Vidal, 17 I. Vila, 11 
R. Vilar, 11 T. Vine, 30 M. Vogel, 36 I. Volobouev«, 28 G. Volpi, 45 F. Wiirthwein, 9 P. Wagner, 44 R.G. Wagner, 2 
R.L. Wagner, 17 J. Wagncr-Kuhr, 26 W. Wagner, 26 T. Wakisaka, 40 R. Wallny, 8 S.M. Wang, 1 A. Warburton, 33 
D. Waters, 30 M. Weinberger, 52 W.C. Wester III, 17 B. Whitehouse, 55 D. Whiteson 6 , 44 A.B. Wicklund, 2 
E. Wicklund, 17 G. Williams, 33 H.H. Williams, 44 P. Wilson, 17 B.L. Winer, 38 P. WitticM, 17 S. Wolbers, 17 
C. Wolfe, 13 T. Wright, 34 X. Wu, 20 S.M. Wynne, 29 A. Yagil, 9 K. Yamamoto, 40 J. Yamaoka, 51 T. Yamashita, 39 
C. Yang, 59 U.K. Yang" 1 , 13 Y.C. Yang, 27 W.M. Yao, 28 G.P. Yeh, 17 J. Yoh, 17 K. Yorita, 13 T. Yoshida, 40 G.B. Yu, 48 
I. Yu, 27 S.S. Yu, 17 J.C. Yun, 17 L. Zancllo, 50 A. Zanctti, 53 I. Zaw, 22 X. Zhang, 24 Y. Zhcng fc , 8 and S. Zucchclli 5 

(CDF Collaboration*) 

1 Institute of Physics, Accidentia Sinica, Taipei, Taiwan 11529, Republic of China 
2 Argonne National Laboratory, Argonne, Illinois 604-39 
"institut de Fisica d'Altes Energies, Universitat Autonoma de Barcelona, E-08193, Bellaterra (Barcelona), Spain 

^Baylor University, Waco, Texas 76798 
5 Istituto Nazionale di Fisica Nucleare, University of Bologna, 1-40127 Bologna, Italy 
''Brandeis University, Waltham, Massachusetts 02254 
7 University of California, Davis, Davis, California 95616 
8 University of California, Los Angeles, Los Angeles, California 90024 
9 University of California, San Diego, La Jolla, California 92093 
10 University of California, Santa Barbara, Santa Barbara, California 93106 
11 Instituto de Fisica de Cantabria, CSIC- University of Cantabria, 39005 Santander, Spain 
12 Carnegie Mellon University, Pittsburgh, PA 15213 
1:1 Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637 
1J> Comenius University, 842 48 Bratislava, Slovakia; Institute of Experimental Physics, 040 01 Kosice, Slovakia 
15 Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 

16 Duke University, Durham, North Carolina 27708 
17 Fermi National Accelerator Laboratory, Batavia, Illinois 60510 
18 University of Florida, Gainesville, Florida 32611 
19 Laboratori Nazionali di Frascati, Istituto Nazionale di Fisica Nucleare, 1-00044 Frascati, Italy 
20 University of Geneva, CH-1211 Geneva 4, Switzerland 
21 Glasgow University, Glasgow G12 8QQ, United Kingdom 

22 Harvard University, Cambridge, Massachusetts 02138 
23 Division of High Energy Physics, Department of Physics, 
University of Helsinki and Helsinki Institute of Physics, FIN-00014, Helsinki, Finland 
2J> University of Illinois, Urbana, Illinois 61801 



3 



25 The Johns Hopkins University, Baltimore, Maryland 21218 
2b Institut fur Experimentelle Kernphysik, Universitdt Karlsruhe, 76128 Karlsruhe, Germany 
27 Center for High Energy Physics: Kyungpook National University, 
Daegu 702-701, Korea; Seoul National University, Seoul 151-742, 
Korea; Sungkyunkwan University, Suwon 440-746, 
Korea; Korea Institute of Science and Technology Information, Daejeon, 
305-806, Korea; Chonnara National University, Gwangju, 500-757, Korea 
28 Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, California 94720 
29 University of Liverpool, Liverpool L69 7ZE, United Kingdom 
30 University College London, London WC1E 6BT, United Kingdom 
31 Centro de Investigaciones Energeticas Medioambientales y Tecnologicas, E-28040 Madrid, Spain 
32 Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 

33 Institute of Particle Physics: McGill University, Montreal, 
Canada H3A 2T8; and University of Toronto, Toronto, Canada M5S 1A7 
4 University of Michigan, Ann Arbor, Michigan 48109 
35 Michigan State University, East Lansing, Michigan 48824 
36 University of New Mexico, Albuquerque, New Mexico 87131 
37 Northwestern University, Evanston, Illinois 60208 
38 The Ohio State University, Columbus, Ohio 43210 
39 Okayama University, Okayama 700-8530, Japan 
J, ° Osaka City University, Osaka 588, Japan 
41 University of Oxford, Oxford 0X1 3RH, United Kingdom 
42 University of Padova, Istituto Nazionale di Fisica Nucleare, 
Sezione di Padova-Trento, 1-35131 Padova, Italy 
43 LPNHE, Umversite Pierre et Mane Curie/IN2P3-CNRS, UMR7585, Pans, F-75252 France 
44 University of Pennsylvania, Philadelphia, Pennsylvania 19104 
4 '' Istituto Nazionale di Fisica Nucleare Pisa, Universities of Pisa, 
Siena and Scuola Normale Superiore, 1-56127 Pisa, Italy 
46 University of Pittsburgh, Pittsburgh, Pennsylvania 15260 
47 Purdue University, West Lafayette, Indiana 47907 
48 University of Rochester, Rochester, New York 14627 
49 The Rockefeller University, New York, New York 10021 
j0 Istituto Nazionale di Fisica Nucleare, Sezione di Roma 1, 
University of Rome "La Sapienza," 1-00185 Roma, Italy 

51 Rutgers University, Piscataway, New Jersey 08855 
52 Texas A&M University, College Station, Texas 77843 
53 Istituto Nazionale di Fisica Nucleare, University of Trieste/ Udine, Italy 
54 University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
55 Tufts University, Medford, Massachusetts 02155 

56 Waseda University, Tokyo 169, Japan 
57 Wayne State University, Detroit, Michigan 48201 
58 University of Wisconsin, Madison, Wisconsin 53706 
59 Yale University, New Haven, Connecticut 06520 
(Dated: March 29, 2008) 

We report on a search for the process pp — > 7 + W/Z with W/Z — > qq in events containing two jets 
and a photon at the center-of-mass energy ^/s = 1.96 TeV, using 184 pb -1 of data collected by the 
CDF II detector. A neural network event selection has been developed to optimize the rejection of 
the large QCD production background; it is shown that this method gives a significant improvement 
in both signal-to-noise ratio and signal sensitivity, as compared with an event selection based on 
conventional cuts. An upper limit is presented for the 7 + W/Z production cross section with the 
W and Z decaying hadronically. 

PACS numbers: 



* Deceased 

tWith visitors from a University of Athens, 15784 Athens, Greece, 
''Chinese Academy of Sciences, Beijing 100864, China, c University 
of Bristol, Bristol BS8 1TL, United Kingdom, d University Libre 
de Bruxelles, B-1050 Brussels, Belgium, e University of California 
Irvine, Irvine, CA 92697, ^University of California Santa Cruz, 



Santa Cruz, CA 95064, 9 Cornell University, Ithaca, NY 14853, 
^University of Cyprus, Nicosia CY-1678, Cyprus, 'University Col- 
lege Dublin, Dublin 4, Ireland, ^University of Edinburgh, Edin- 
burgh EH9 3JZ, United Kingdom, ^University of Heidelberg, D- 
69120 Heidelberg, Germany, ; Universidad Iberoamericana, Mexico 



4 



I. INTRODUCTION 

The identification of gauge boson hadronic decays is 
extremely challenging at hadron colliders, since a small 
two-jet resonance needs to be extracted from a huge 
QCD multi-jet background. At the Tevatron only the 
favorable circumstance of W's generated in top quark 
decays has allowed for a successful identification of the 
W hadronic resonance 1\. Nevertheless, the ability to 
extract hadronic resonances submerged in a large QCD 
background is of paramount importance in the search 
of new particles with dominantly hadronic decays. The 
most important example is the Higgs boson for which no 
direct evidence has yet been observed. 

At the Tevatron, one of the most promising signatures 
for the Higgs observation is the associated production 
with a W(Z), where the Higgs decays into two jets 0. 
However, at the center of mass energy \fs = 1.96 TeV, 
the standard model (SM) Higgs boson cross section is 
much smaller than that for the non-resonant W +jj pro- 
duction, and thus, sophisticated techniques are needed to 
suppress the QCD background while maintaining a high 
signal detection efficiency. 

In this respect, identification of dijet resonances of the 
W and Z bosons provides an important test bench for 
developing such techniques, due to the high statistical 
sample that can be collected and the fact that their char- 
acteristics are well known. In addition, a highly popu- 
lated W/Z boson dijet mass peak is an excellent tool to 
constrain the jet energy scale and also to improve the di- 
jet mass resolution, two essential ingredients for precision 
measurement of signatures with jets in the final state. 

At hadron colliders, a mass peak from W(Z) — > jj was 
reconstructed in the inclusive dijet events by the UA2 
collaboration [3[ at </s = 630 GeV. With a signal over 
background ratio (S/B) of about 1/35, about 5000 events 
were observed. At \/s — 1.96 TeV, the QCD dijet pro- 
duction cross section increases by approximately a factor 
35 for 20 GeV jets, making the production rate too high 
to be handled by the data acquisition system. However, 
this is not the case when the W(Z) is produced in asso- 
ciation with another gauge boson (j,W,Z). 

Because the 7 + W(Z) cross section is one order of 
magnitude higher than the heavy diboson production 
WW + WZ, these events offer in principle the best op- 
portunity to identify the W{Z) — > jj resonance. 

In addition, the diboson production with a photon is 
interesting in its own right. In fact, the 7 + W(Z) pro- 
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duction is directly correlated to the non-Abelian charac- 
ter of the electroweak theory, and is sensitive to physics 
beyond the standard model through enhancement of the 
trilincar WW"f coupling and possible contributions of the 
ZZj and Z77 couplings forbidden in the standard model. 
Although such effects have already been searched for in 
the leptonic channels of W(Z)j events the successful 
identification of such events also in the hadronic channels 
could concur for an even more stringent test of the SM 
in this sector. 



A. Analysis Overview 

In this paper we report on a search for W(Z) decaying 
into two jets based on a sample of 7 + jj data collected 
with the CDF II detector between July 2003 and Septem- 
ber 2004 H, corresponding to an integrated luminosity 
of 184±7 pb _1 j|. In a previous study of this signature 
performed by the CDF collaboration at \/s = 1.8 GeV 
and using 90 pb" 1 of data Q, a significance (S/y/S + B) 
of 0.3 was achieved, with a S/B of about 1/100. In the 
study reported here, in addition to an improved online 
event selection, a neural network based technique is em- 
ployed to enhance the significance. 

The expected shape of the WiZ) mass distribution 

(jrijj^ Z ) is derived from simulated SM signal events. The 
shape of the background is determined directly from the 
data by fitting the observed dijet mass distribution (rrijj) 
in the control region, i.e. excluding the part of the rrijj 
spectrum around the W/Z boson mass value where the 
signal is expected to be visible (signal region). 

Because of the steeply falling behavior of the rrijj dis- 
tribution, it is important to have unbiased control regions 
both below and above the signal region to obtain an ac- 
curate description of the background. Extreme care is 
taken in choosing the online and offline selection cuts in 
order to not deplete the control region at low values of 
rrijj. In fact, since such region has the biggest weight in 
the fit, it ultimately determines the accuracy with which 
the background estimate can be determined. Such accu- 
racy is particularly crucial in cases with very low S/B 
ratio, like the search reported in this article. The excess 
in the signal region over the smooth background - if con- 
sistent with the SM signal shape - can then be attributed 
to W(Z) decaying into jets. 

This paper is organized as follows. In section HH a 
description of the processes involved in the W(Z)j pro- 
duction is provided as well as the SM cross section pre- 
dictions. Detector and trigger descriptions follow in Sec- 
tions IIIII and IIVI In Sections [V] and I VII event selection 
criteria and expected event yield, along with their sys- 
tematic uncertainties, are outlined. The neural network 
based selection and its performance is described in Sec- 
tion IVII1 Sections IVIIII and IIXI discuss the results fol- 
lowed by the conclusions. 
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II. STANDARD MODEL PREDICTION FOR 
THE W(Z)j CROSS SECTION 

The tree-level Feynman diagrams for W7 and Zj pro- 
duction are shown in Fig. 1. Figure 1(a) and Fig. 1(b) 
show the ^-channel and u-channel W(Z) production re- 
spectively, where a photon is radiated from one of the 
incoming quarks. Figure 1(c) and Fig. 1(d) show the pro- 
cesses where a photon is radiated from the decay quarks 
of the W(Z) boson. In these latter cases the W(Z) boson 
resonance cannot be reconstructed from the two-body 
mass of the final quarks. The final state of these pro- 
cesses is very similar from both kinematic and topologi- 
cal standpoints to some components of the background in 
our sample (Sec. IVIIA]) . Because our analysis cuts have 
a high background rejection power, 1(c) and 1(d) radia- 
tive decays contributions are strongly suppressed in the 
sample selected. Finally the process involving the three 
vector boson coupling WWj is shown in Fig. 1(e). 

The pp — > W{Z)"f — > qq'f predictions were determined 
using the PYTHIA Q Monte Carlo (MC) generator. It 
calculates the matrix elements at leading order (LO) and 
in the narrow boson-width approximation in which ra- 
diative boson decay diagrams (Fig. 1(c) and 1(d)) are 
ignored. The QCD initial/final state radiation as well 
as subsequent parton fragmentation and hadronization 
were also provided by PYTHIA . The factorization scale 
Q was set equal to the center of mass energy of the in- 
coming quarks y/i. The CTEQ5L Q parton distribution 
functions were used. 

The PYTHIA calculations were compared to the pre- 
dictions obtained with MADGRAPH [ll(, a tree level ma- 
trix element calculator that, in contrast to pythia, does 
not use the narrow boson-width approximation. The ra- 
diative contribution was suppressed by requiring the in- 
variant mass of the W(Z) di-quark decays to be greater 
than 74(85) GeV. In addition, to avoid collinear emis- 
sion divergences the distance between the photon and the 
quarks in the r]-(j) [lCj space was required to be greater 
than 0.4. The final state observables from the MC sim- 
ulation were also compared. The px and pseudorapid- 
ity distribution of the two outgoing partons and of the 
photon as well as the pair-wise separation, defined as 

A% = V(&-h) 2 + (m-Vj) 2 i = 1,2,3, i + j, 
were in excellent agreement. The resulting madgraph 
W'y (Zj) cross section is 11% higher than the value pre- 
dicted by pythia . The PYTHIA cross section prediction 
is scaled for O(as) QCD contributions (fc-factor) com- 
ing from subprocesses with either virtual gluon loops or 
gluon/quark emissions in the initial state. The magni- 
tude of this correction, averaged over the photon spec- 
trum in the region p T > 10 GeV, is 1.55 for the Wj 
process [l2j and 1.44 for the Zj [l3j]. Including this fc- 
factor the SM prediction for the signal in the kinematic 
region pip > 10 GeV and |7y 7 | < 1.2 is 

°W(Z)~t = V{PP -> Wl) x <B(W -> qq) 

+ a{pp -> Zi) x <B(Z -> qq) = 20.5 ± 2.5 pb. 



The 12% uncertainty accounts for the discrepancy be- 
tween PYTHIA and MADGRAPH cross sections (11%), for 
the fc-factor (3%), the factorization scale (1.5%), and the 
parton distribution function (4.8%) uncertainties. It is 
interesting to notice that in contrast to the inclusive pro- 
duction, where the W cross section is about three times 
larger than the Z cross section Q, for the production of 
the W and Z in association with a photon, the SM pre- 
dicts similar cross sections (awy = 9.9 pb and az-y = 10.6 
pb). 



III. DETECTOR DESCRIPTION 

A detailed description of the CDF II detector can be 
found elsewhere [14| . Here we briefly describe the aspects 
of the detector relevant for this analysis. The track- 
ing system is a magnetic spectrometer consisting of a 
90-cm long cylindrical silicon micro-strip detector sur- 
rounded by a 3.1 m long drift chamber, both immersed 
in a 1.4 T magnetic field. The calorimeter consists of 
an electromagnetic (em) and a hadronic (had) compart- 
ments covering both central (|?7| < 1.1) and forward 
(1.1 < < 3.6) regions. Both calorimeters are seg- 
mented into projective towers. The tower size in the cen- 
tral calorimeter is approximately 0. 11(77) x 15°(0), and 
the resolution is about 13.5%/V-Et © 2% for electrons 
(where Et = Esin8 and E is measured in GeV). Em- 
bedded in the central calorimeter is a a multiwire pro- 
portional chamber (ces), located at a depth of approx- 
imately six radiation lengths where the density of the 
energy deposited by an EM shower is at a maximum. 
Cathode strips and anode wires, with a channel spacing 
between 1.5 and 2 cm, running along the azimuthal and 
the beam line direction respectively provide precise in- 
formation on the electromagnetic shower centroid as well 
as the shower profile in the transverse direction. Another 
wire chamber (cpr) is located between the magnet coil 
and the central calorimeter modules. It measures the 
signals from early showers of electromagnetic particles 
occurring in the coil. The CES and CPR systems are used 
to discriminate prompt photon from multi-photon decay 
products of neutral mesons, 7r°'s, 77's or Ks's. 

The data were collected with a three level trigger sys- 
tem. At level 1 (LI), a simple selection can be made 
based on the presence of tracks above a fixed pr thresh- 
old, on the total energy deposited in the calorimeter, 
or on single calorimeter trigger tower energies (a trig- 
ger tower consists of two calorimeter towers adjacent in 
the z direction). At level 2 (L2), custom built hard- 
ware is used to reconstruct calorimeter energy clusters, 
apply isolation requirements for photons and electrons, 
identify muons, and measure track displacements from 
the primary vertex. At level 3 (L3), events are fully re- 
constructed with the same algorithms used in the offline 
analysis. The transverse energies however are calculated 
using the nominal interaction point, instead of the actual 
event vertex position. 
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. . . x (d) Inner Bremsstrahlung 

(ej s — channel (trihnear coupling) 

FIG. 1: (a,b,e) Feynman diagrams at the tree level for the process qq — » W(Z) r y — > qqy. The s-channel for the Z'y production 
is forbidden in the standard model. (c,d) Radiative W(Z) decays diagrams. A bremsstrahhmg photon is emitted by one of the 
two quarks from the W(Z) decay. 



IV. TRIGGER SELECTION 

In the analysis reported in this article the main source 
of background is the non-resonant QCD 7 + jj produc- 
tion. In addition, a large contribution from three-jet pro- 
duction is also expected. Both of these background pro- 
cesses have rather large event rates. As a consequence, an 
elaborate triggering scheme is needed to reduce their rate 
to levels that can be handled by the current data acquisi- 
tion hardware. The main challenge is to keep the photon 
Pt threshold low enough in order not to bias significantly 
the data rrijj distribution below the signal region. Only 
with this requirement can an accurate determination of 
the background shape be successfully carried out (see 
Sec. II A|) . However, an inclusive photon trigger with a 
low p T threshold results in an unacceptably large rate. 
We designed a trigger taking into account the above con- 
straints. Details of the trigger specifications are outlined 
in the following sections. 



A. Level 1 and level 2 Selection 

At level 1 events with a trigger tower with Et > 8 GeV 
and at least 89% of its energy deposited in the EM sec- 
tion are selected. At level 2, electromagnetic clusters are 
reconstructed combining towers with Et > 7.5 GeV ad- 
jacent to a seed tower. A seed tower must have an Et > 8 
GeV with 89% of its energy deposited in the em calorime- 
ter. Only em clusters with Et > 12 GeV and isolated 
from other deposits of energy are selected. The isolation 
requirement proceeds as follows. The sum of the trans- 
verse energies is determined in a) 8 towers surrounding 
the seed tower and b) all four combinations of ten towers 



in a 4 x 3 region surrounding the seed and one adjacent 
tower. The lowest of these five sums is required to be less 
than 1 GeV. Such a strict isolation requirement provides 
significant rejection against the high-rate neutral meson 
multi-photon-decay background and against photons ra- 
diated by quarks or gluons. To further reduce the L2 
output rate the presence of a significant hadronic activ- 
ity were added on top of the photon requirement. The 
L2 hardware jet finder was exploited to identify clusters 
of energetic towers where a nearest neighbor algorithm 
with a seed tower threshold of 3 GeV is used. The trig- 
ger requires the presence of at least two such L2 clusters, 
one of which corresponds to the photon, with the seed in 
the region \rj\ < 1.78. To maximize efficiency for low Et 
jets we apply no explicit requirement on cluster energy. 
Instead, the total transverse energy of the calorimeter 
trigger towers TiEt, excluding the photon candidate en- 
ergy, is required to be greater than 20 GeV. The trigger 
rate reduction brought about by these extra cuts allows 
the photon Et threshold to be set as low as 12 GeV. 



B. Level 3 Selection 

At level 3, em clusters are formed by combining towers 
with more than 2 GeV of energy with their two nearest 
neighbors in pseudorapidity. Only clusters with 95% of 
their energy in the electromagnetic calorimeter are se- 
lected. Positions and transverse profiles of em cluster 
showers are determined using the CES detector. Eleven 
strips (wires) around the most energetic strips (wires) 
are grouped to form a CES cluster. To avoid spurious 
clusters made up by noisy channels, at least two strips 
(wires) in each cluster are required to be above thresh- 
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old. This solution is more efficient than just requiring one 
strip (wire) with high energy, as was done in Run I [la ]. 
The precise position of an em cluster is determined us- 
ing the centroid of the most energetic CES cluster inside 
the EM cluster towers. The position resolution achieved 
using this method is about 2 mm for a single particle 
shower. The CES cluster centroid is also required not 
to be close to the edges of the ces where the detector is 
not fully efficient. In particular it has to be within 21 
cm from the center of the tower in azimuthal direction 
(Xces) and within 9 cm < \z\ < 217 cm along the beam 
direction (Zqes)- A calorimeter cluster isolation energy 
£™ is defined at this level as the total transverse energy 
inside a cone of radius R=0.4 in ij — <p space, centered 
at the CES cluster position, but excluding the cluster 
energy. A cut of E^ ° < 1 GeV is applied. The profile 
of the cluster is compared with a single em particle pro- 
file as measured in test beam and \ 2s quantifying the 
"similarity" are formed in both the azimuthal and longi- 
tudinal directions [lj|. The average of these two x 2 s are 
required to be less than 20. No explicit requirements on 
jets are implemented at L3. A summary of the trigger 
cuts is reported in Table. 1. 



TABLE I: Summary of the requirements implemented in the 
trigger at different levels. At level 1 and level 3 only cuts 
on photon related quantities are implemented. At level 2 
requirements on hadronic clusters are present as well. 

Level 1 

Trigger Tower E T > 8 GeV 

Trigger Tower E HAD /E EM < 0-125 

Level 2 - Photon Cuts 
L2 em Cluster E T > 12 GeV 

L2 em Cluster Eh ad / Eem < 0.125 
L2 em Cluster \r]\ < 1.2 

L2 em Cluster E™° < 1.0 GeV 



L2 J2 e t 
L2 Jet 
L2 Jet |r?| 



Level 2 - Jet Cuts 

> 20 + Vt GeV 
> 1 
< 1.78 



Level 3 

L3 EM Cluster E T 
L3 EM Cluster E had /E E m 
L3 em Cluster E^ so 
L3 em Cluster Xces 
L3 EM Cluster \X C es\ 
L3 em Cluster \Zces\ 



> 12 GeV 

< 0.05 
< 1.0 GeV 

< 20 
< 21 cm 
9 < z < 217 cm 



V. EVENT SELECTION 

The events selected online are processed offline tak- 
ing into account the updated calorimeter calibration, the 
tracker alignment constants, and the measured beam po- 
sition in the data. The primary vertex location is de- 
termined by iteratively fitting the tracks to a common 



point. In case more than one vertex is reconstructed due 
to multiple pp interactions in the same bunch crossing, 
the primary vertex of the event is considered that whose 
associated tracks have the highest sum of transverse en- 
ergy. The transverse energies are then determined with 
respect to this interaction. In the following the offline 
event selection is described. 



A. Photon Selection 

To eliminate the cosmic ray contamination from the 
sample, the total missing transverse energy [l0( is re- 
quired to be less than 80% of the transverse energy of 
the photon candidate. The primary event vertex posi- 
tion along the beam direction is required to be within 60 
cm from the center of the detector. Only events with an 
EM cluster with Et > 12 GeV are selected. The cluster 
position determined in the CES detector is restricted to 
\X C es\ < 17 cm and 14 cm < \Z C es\ < 217 cm. These 
fiducial cuts ensure the EM shower is contained inside 
the CES detector boundaries, allowing an accurate re- 
construction of its transverse profiles. The isolation cut 
applied at the trigger level is refined offline where the 
transverse energy in a cone R = 0.4 around the EM clus- 
ter, calculated using the event vertex, is required to be 
less than 1 GeV excluding the photon transverse energy. 
The photon energy is corrected in average for the con- 
tributions of multiple pp interactions in the same bunch 
crossing (pile-up events) and for the photon em shower 
leakage into neighboring towers. In addition, the isola- 
tion requirement is reinforced by rejecting photon candi- 
dates with a reconstructed track pointing to it. Photons 
converted into e + e~ pairs in the tracking volume or in 
the beam pipe, about 14% of all photons emerging from 
the interaction point, are also rejected by this cut. The 
CES shower shape is compared to the one generated by a 
single em particle profile with the same technique used 
at L3. A similar Xces < ^0 cut is thus applied. Photon 
candidates with a second CES cluster inside the associ- 
ated EM cluster and with energy above 1 GeV are also 
rejected to suppress the multi-photon background. The 
efficiencies of these cuts in selecting prompt photons are 
described in Section IVIBI 



B. Photon Background Subtraction 

The photon candidates passing the above requirements 
are still contaminated by multi-photons from neutral me- 
son decay. Two independent techniques are employed 
to subtract this multi-photon background on a statis- 
tical basis. The first one ("profile method") exploits 
the difference in Xces °f the ^ w0 components. Low px 
prompt photons are expected to have a smaller Xces 
than multi-photons which have a broader em shower pro- 
file. However, this method is not useful for EM clusters 
with pt > 35 GeV: at such energies multi-photons are 
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too collimated to produce electromagnetic showers that 
are dctcctably broader than single photon. The second 
technique ("conversion method") [1 51 ] exploits instead 
the different conversion probability of single and multi- 
ple photons when they pass through the magnet coil, and 
it is approximately independent of pr- Such conversions 
are detected in the CPR detector. For both methods the 
prompt photon content of the sample is given by: 



N, 



total i 



where e 7 and ef, are respectively the efficiencies for 
prompt and multiple photons to pass a fixed Xces cu ^ 
(CPR pulse height cut) in the case of the profile (con- 
version) method. Such efficiencies are determined using 
both real data and simulated control samples as detailed 
in [161 ] . The number of photon candidates in the sample 
is N total and e is the fraction of these candidates passing 
the cuts. The two methods provide a consistent estimate 
of the prompt photon content. In the following, for pho- 
ton background subtraction, we determine the prompt 
photon content using the profile (conversion) method to 
photon candidates with px < 35 GeV {pr > 35 GeV ). 
All the event distributions, including rrijj, are accord- 
ingly weighted to subtract the multiphoton background. 
The ratio of the number of prompt photons to the num- 
ber of photon candidates in the sample after the event 
selection is shown in Fig.H 




10 15 20 25 30 35 40 45 

Photon Candidate p (GeV/c) 



FIG. 2: Prompt photon fraction (number of prompt pho- 
ton/number of candidate photons) in the data sample as a 
function of the photon candidate pr- Only statistical errors 
are reported. 



C. Jet Selection 

Hadronic jets are identified using an iterative cone clus- 
tering algorithm [13] with a cone radius R=0.7. Based 



on simulations of jet fragmentation and of calorimeter 
response to hadrons the raw Et of the jets are corrected 
for [3: 

(i) the non-linear and non-uniform response of the 
calorimeter; 

(ii) the undetected energy falling into uninstrumented 
regions of the detector; 

(iii) the energy coming in average from multiple pp in- 
teractions occurring in the same bunch crossing and 
the underlying event contribution; 

(iv) the energy of low momentum charged particles that 
do not reach the calorimeter; 

(v) the average energy loss due to particles falling out- 
side the jet-clustering cone. 

The jet corrections depend on the px of the jet, its pseu- 
dorapidity, and on the number of vertices in the event. 
They amount, on average, to 25%(15%) of the jet energy 
for 15(50) GeV jets. In this analysis only events with 
two jets of Et > 15 GeV and containing no additional 
jet with Et > 10 GeV are selected. The additional jet 
veto is introduced both to reduce the QCD background 
and to improve the W/Z dijet mass resolution by remov- 
ing W{Z)^i events with hard gluon radiation. 



VI. SELECTION EFFICIENCY AND SIGNAL 
YIELD 

In this section the trigger and offline requirement ef- 
ficiencies in selecting ^(W/Z) — > jqq signal events are 
calculated. The trigger efficiency is calculated for events 
satisfying all the offline selection criteria. In turn, the 
offline selection efficiency is evaluated using simulated 
W(Z)j events. 



A. Trigger Efficiency 



It is convenient to break up the trigger efficiency 
in two components: 1) the photon selection efficiency e 7 



t 

1 rg 



and 2) the efficiency related to hadronic cluster require- 
ments e{f q s (see Tab. HVB]) . 



The . 



trg 
.7 
trg 



value is calculated as follows. First, it is eval- 



uated relative to a control sample collected by a trigger 
with looser photon cuts (including a lower pt threshold) . 
Then the efficiency of this control sample is measured 
using a sample of "unbiased" photon candidates, i.e. a 
sample where they have not been used to trigger the data 
set. The product of these two contributions is shown in 
Fig. [3J this gives the photon candidate trigger efficiency. 
The value at the plateau reflects the online/offline isola- 
tion energy differences while the low pt turn-on is deter- 
mined by the trigger threshold energy smearing. The fi- 
nal prompt photon trigger efficiency ej rg is determined by 
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applying the photon background subtraction described in 
Sec. IV Bl to the plot in Fig. [3J For prompt photons the 
plateau level increases to 85% as they are more likely to 
pass the isolation cuts than the multi-photons. 



.2 1 
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FIG. 3: Trigger efficiency for photon candidates as a function 
of pt ■ 



The 



Jets 
-trg 



is evaluated using the MC signal sample 



described in Section [TTJ The simulation of the CDF 
calorimeter has been tuned to reproduce the response 
measured in collider data [HI]. The energies in the trig- 
ger towers, the L2 Et> an d the L2 jets quantities were 
estimated using the online algorithms applied to the of- 
fline calorimeter tower energies. The accuracy of such 
estimates was checked in real data against the actual on- 
line measurements, and for the quantities used in the 
online selection the agreement was found to be within 
1%. The fraction of MC signal events passing the L2 jet 
requirements is e^g = 0.93. 

In conclusion, the combined trigger efficiency in select- 
ing j(W/Z) —>■ 'jqq event is 



e tra 



c 7 



Jets 



trg trg 



0.76 ±0.01. 



B. Acceptance and Selection Efficiency 



The acceptance and efficiency of the offline event se- 
lection is estimated by applying sequentially the cuts de- 
scribed in Section [ill to the MC signal sample. In Ta- 
ble pi] the offline cut relative efficiencies, defined as the 
fraction of events passing a cut after having passed all 
the previous cuts, are reported. The MC simulation ac- 
ceptance - the fraction of generated events containing an 
EM cluster of Et > 12 GeV - reflects the choice of the 
Pt photon generation cut (10 GeV). A lower cut at the 
generation level is needed to avoid threshold bias brought 
about by the finite detector resolution. The photon ge- 
ometric acceptance includes the pseudorapidity selection 



as well as the Xq es and Zqes cuts. The accuracy of the 
efficiencies reported in Table [TT] depends upon the preci- 
sion of the detector simulation in reproducing the data. 
The electromagnetic particle response in the simulation is 
checked using electrons from Z — > ee and W — > ev decays 
(a large sample of pure prompt photons is not available 
in the data). This comparison is used to estimate the 
systematic uncertainties of the selection efficiencies. An 
account of these studies is given next: 

(i) Z vertex Cut: The shape of the luminous region in 
real data was determined by fitting the vertex po- 
sition in minimum bias events. The signal vertex 
position is simulated according to this distribution. 
The fraction of events within \z\ < 60 cm in MC 
simulation matches the data within 0.5%. 

(ii) Missing Energy Cut: A change of the E™ ss /p T by 
10% resulted in a 2% change in the selection effi- 
ciency, which is assigned as systematic uncertainty. 

(iii) had/em Ratio: The fraction of the unbiased elec- 
tron from Z-boson decays which pass the had/em 
cut in simulated and data events agrees within 1%. 
We assume the same difference holds for photons, 
whose shower starts deeper in the calorimeter, and 
assigning a 1% systematic uncertainty. 



(iv) Calorimeter Isolation: The amount of energy sur- 
rounding a prompt photon em deposition is deter- 
mined by em shower leakage outside the cluster and 
by underlying and multiple interaction events. The 
accuracy of the simulation of the isolation cut mea- 
surement has been evaluated using cones of R = 0.4 
randomly placed in the photon fiducial region. The 
energy collected in these cones can be considered 
an approximation of the isolation energy measured 
around em clusters. The fraction of such cones pass- 
ing the isolation cut (Et < 1 GeV) in simulated 
W — > ev events was found 3±2% higher than in 
the data. A correction factor 0.97 is applied to the 
MC isolation efficiency to account for the observed 

(1) discrepancy. The 2% uncertainty is included in the 

systematic errors. 

(v) Track Isolation: The track isolation efficiencies in 
simulated and real data events were found to be 
consistent within 2%. The photon conversion 7 — > 
e~e~ rate is used to tune the detector simulation 
for the amount of material present in front of the 
calorimeter. The uncertainty of the track isolation 
efficiency includes any remaining deficiency in the 
material simulation. 

(vi) CES x 2 ■' Photon and electron EM shower profile are 
simulated using the information collected during the 
single electron test beam. As a consequence, for the 
efficiency of the CES x 2 cut, a very good agreement 
(within 0.2%) is observed between simulated and 
real data Z — > ee events. 
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(vii) CES Cluster Isolation: The CES cluster activity 
around Z-boson decay electrons in MC simulation 
was found to match the data at the level of 3%. 

(viii) Jet Cuts: To assess the uncertainty on the jet cut 
efficiency the jet energy scale of all jets is shifted 
by one standard deviation (about 8(4)% for jets of 
E T =15(50) GeV [3). This results in a 7% relative 
change on the selection efficiency that we set as sys- 
tematic uncertainty. This is the dominant source of 
systematics. 



TABLE II: Summary of the event selection cuts and their 
relative efficiency. 



Analysis Cuts 


Efficiency (%) 


MC Simulation Acceptance 


62.1±0.1 


Photon Geometric Acceptance 


60.7±1.0 


\z v tx\ < 60 cm 


96.1±0.5 


Missing E T Cut 


90.6 ±2.0 


Total Acceptance: Akin = 0.33 ± 0.01 


had/em Ratio 


94.5±1.0 


Calorimeter Isolation 


80.8±2.0 


Track Isolation 


80.2 ±2.0 


CES x 2 


99.6±0.2 


CES Cluster Isolation 


94.8±3.0 


Total Photon Identification Efficien 


cy: cp h = 0.58 ± 0.05 


Jet 3 p T < 10 GeV Cut 


54.6 


Jet 2 p T > 15 GeV Cut 


82.8 


Total Jet Selection Efficiency: 6j et 


= 0.45 ±0.03 



Combining all the contributions in Table |TT] (Akin , £ph , 
and ej e t) with the trigger efficiency et rg , our estimate of 
the total signal selection efficiency is 

e = e trg ■ A kln ■ e ph ■ ejet = 0.065 ± 0.006. (2) 

C. Signal Dijet Mass Distribution and Signal Yield 

The mass distribution of the two leading jet system 
for the simulated signal events passing all the Table |TT] 
selection criteria is reported in Fig. [4] along with the in- 
dividual jW and jZ contributions. Both the W and Z 
mass distributions have non-Gaussian tails arising from 
initial and final state gluon radiation. For the Z we no- 
tice a larger low mass tail due to the higher - on average - 
quark momenta compared to the W quarks. In the range 
between 60 and 120 GeV the signal can be adequately de- 
scribed by a single Gaussian with a mean value of 87.2 
GeV and a width of 12.5 GeV. This shape is used to 
extract the signal from the data. The dijet mass resolu- 
tion (AM/M), estimated by fitting a Gaussian function 
around the W and Z peaks, is 12% for both gauge bosons. 
This is consistent with other MC dijet mass resolution 
studies [|| . The expected number of signal events in the 
sample is given by N = e x a^w/z x A where e = 0.065 is 



W/Z Peak (n=87.2 o=12.5) 




Dijet Mass (GeV) 



FIG. 4: Normalized dijet invariant mass (solid line) distribu- 
tion of the two leading jets of the selected *y(W/Z) — > 'yqq 
MC events. The individual contributions from the W and Z 
bosons (dashed lines) are shown. The fit results are reported 
in the inset. 

the selection efficiency (without any mass window cuts), 
v-yW/z = 20.5 pb is the SM cross section (reported in 
Section [IT}, and C = 184 pb -1 is the total integrated 
luminosity of the sample. In the dijet mass window 60 
< rrijj < 120 GeV, 227 signal events are expected among 
the 42462 events present in the data. This corresponds 
to a signal over background ratio (S/B) of 1/187. For 
the current data set, the statistical significance - defined 
as S/VS + B - is 1.1. 

In the following we show how the use of a neural net- 
work in the selection process can substantially enhance 
the sensitivity of the analysis. 

VII. ADVANCED EVENT SELECTION 

The basic event selection described in Sec. [V] has a 
rather mild discrimination power. However, the 15 GeV 
jet px threshold cannot be increased since it would de- 
plete the low end of the dijet mass spectrum. Similarly, 
the rejection of the extra jet activity is meant more to 
improve the dijet mass resolution rather than suppressing 
the background. Nevertheless, the kinematic and topo- 
logical distributions of the final state in signal and back- 
ground events exhibit some differences that can poten- 
tially be useful in enhancing the sensitivity of this anal- 
ysis. In fact: 

(i) In signal events the W(Z) boson has a low pr (since 
p^^ ~ p^,). As a consequence the two jets are ba- 
sically back-to-back with approximately the same 
energy, and the jet and photon directions are not 
correlated. In contrast the dominant background 
(7 + jj events) comes either from a qg — > qj pro- 
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duction, where the quark balancing the photon ra- 
diates a gluon, or from a qq/qg production, where 
one of the two outgoing quarks radiates a photon. 
In both cases, the radiated gluon/photon tends to 
be collinear with the radiating quark. Thus, the 
photon is either along or in opposite direction to 
the leading jet in the event. 

(ii) In 7 + jj events the two leading jets are typically 
a quark and a gluon jet. This is also true for dijet 
production which is dominated at low px by quark- 
gluon scattering. For signal events instead the two 
leading jets are both quark jets. 

(iii) The signal is characterized by the production of two 
colorless gauge bosons that constrains the initial 
and final state in a particular (color singlet) con- 
figuration. The QCD background involves instead 
quarks and gluons with multiple color connections 
resulting in higher color radiation. 

Hence, it is clear that the production of the signal and the 
background events differs in many ways. However, it was 
not possible to identify a set of selection criteria able to 
adequately discriminate between signal and background 
while keeping an acceptable signal yield. This is shown 
in Fig. [5] where the signal and background distributions 
for a few observables are compared. For these reasons we 
developed an artificial neural network (ANN) selection to 
exploit subtle differences and variable correlations. The 
ANN selection is applied to the events that have already 
passed the simple kinematic cuts described in Table HU 
The structure of the ANN along with its performance is 
described next. 



A. Neural Network Selection 

In this analysis we employed the jetnet fl9| software 
package to construct a feed-forward network [20]. The 
architecture of the network consists of one intermediate 
(hidden) layer and a single output node. For the net- 
work output Nout a target value of 1 for the signal 
and for the background is chosen. The training for 
the signal recognition is performed using as a template 
"/(W/Z) — > g<?7 events generated by PYTHIA (SeaQl]). As 
background template instead, a subsample of real data 
events is used. In fact it is not trivial to simulate properly 
the QCD 7 + jj production because of the interplay be- 
tween the components associated to the hard process (de- 
termined by matrix element calculations) and the compo- 
nents generated by the development of the hard partons 
(described by parton shower calculations) [21| . In addi- 
tion, further complications arise from NLO effects that 
cannot be neglected for an accurate determination of the 
shape of the observable distributions [22j , a key ingredi- 
ent in an ANN training. Considering that less than 0.6% 
of the data are signal events - based on predicted pro- 
duction rates - data provide an excellent approximation 



for background distributions. Only data events in the 
60 < rrijj <120 GeV signal mass window are considered 
in the ANN training. 



B. Variable Selection and Neural Network Training 
Tuning 

We consider a set of 19 input variables (or nodes) re- 
lated to the signal and background differences outlined 
above. The selected variables emphasize event and jet 
topologies, rather than absolute kinematic values of the 
final state objects. This is done to preserve as much as 
possible the shape of the rrijj spectrum. The list of the 
ANN input nodes are given in Table IIIII along with their 
definitions. In order to improve the performance of an 
ANN, it is usually advisable to remove fully correlated 
variables from the set of input nodes. To identify among 
our 19 variables the redundant ones, we develop a "rank- 
ing" method that proceeds as follows. 

First, the most discriminating variable is determined 
by comparing the performance of 19 ANN's having each 
variable in Table Hill as a single input node. 

The ratio S/\/S + B for a signal acceptance of 75% is 
used as a figure of merit. 

Second, two input node ANN's are built. They have 
as a first input node the variable found before and as 
a second node one of the remaining variables. The sec- 
ond best variable, defined as the property that provides 
the best discrimination power when paired with the first 
variable, is determined by comparing the significance of 
these ANN's. 

The procedure is repeated, determining at each step 
the variable which, in conjunction with the best set of 
variables found in the previous step, forms the best per- 
forming ANN. At the end, when all the variables are 
considered, an ordered list of properties is generated. In 
Table IIIII the properties are listed in the order resulting 
from this procedure. The highest and lowest significance 
of the ANN's built at the step k (k = 1, .., 19) is shown 
in Fig. [5] The ANN's discriminating power improves 
with the number of input nodes until the properties that 
are subsequently added become strongly correlated with 
those already considered. At this point a plateau in per- 
formance is reached. In our case such a plateau appears 
at about k — 10. Hence, only the first ten properties 
listed in Table IIIII are used as input nodes in the final 
ANN. 

The number of nodes, Nh, in the hidden layer is set 
to 17. Several ANN's with Nh from 11 to 30 were com- 
pared and no significant differences in performance were 
observed. 
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Property 


Description 




r\ separation between the two leading jets 


IL trk 


Number of tracks inside a cone of size 0.5 


in rj — <f> around the two leading jets 


Mj 2 1 E j2 


Mass over energy ratio of the second jet (M = \/E 2 — P 2 ) 


jets 
Imax 


Maximum rj of the two leading jets 



"Intrajet Energy" defined as f2 = [Y^Et 



where ^ Et is the Et scalar sum of the calorimeter towers 
in the pseudorapidity region (r/ DW — 0.3) < 77 < (r/ UP + 0.3) 
and AL = rf p - r, DW + 0.6 with if w = min(^' etl , rf et2 , rf) 

and 77 = max(i7 , ' eil , rf et2 1 rf). 
The energies of the photon and the two jets are uncorrected 



E 



•jet l 



El 



-JiTT 
T 



E; 



JeTT 



j j 



Azimuthal angle between the two jets 



max A$ 



Maximum azimuthal separation between photon and jets 



min A$, 



Minimum azimuthal separation between photon and jets 



Sphericity 



S = 3/2 • (Q 2 + Q 3 ) with < S < 1 



min Arjj~ 



Minimum 77 separation between photon and jets 



max Arjj 



Maximum rj separation between photon and jets 



■yW 



Azimuthal separation between the photon and the jetl-jet2 system 



V32 



Pseudo-rapidity of the second jet 

JeTT 



AE 3 T 3 



Transverse energy difference between jets 



(3 of the jetl-jet2 system 



Aplanarity 



A = 3/2 ■ Q 3 with < A < 0.5 



cos 8* cosine of the angle 6* between the photon and the leading jet 

directions calculated in the 7-jet reference frame 



Ar; 7W / 



77 separation between the photon and the jet-jet system 



TABLE III: Definition of the properties considered as input nodes for the neural network. The sphericity and aplanarity are 
defined after [H. 



Neural Network Output and Improvement in 
Significance 



D. Dijet mass spectrum 



After the training, the ANN can be seen as a function 
associating a real number 0.0 < Nout < 1-0 to each 
event. The Nout distributions for the signal and back- 
ground samples are shown in Fig. [7] Selecting events 
above some ANN output value Nqut clearly enhances 
the signal sensitivity of the sample. 

In Fig. [5] the signal (background) efficiency is shown 
as a function of Ncut- In order not to deplete the sig- 
nal yield too much we set Ncut = 0.6. For this value 
the signal efficiency of the ANN selection, £nn, is 72%. 
After the ANN selection the expected number of signal 
events is S — 164, while 11691 data events remain in the 
60 < Mjj < 120 GeV mass window. This corresponds to 
an S/B — 1/71 with a significance S/\J S + B = 1.51, an 
improvement of 163% (37%) in S/B (S/VS + B) o ver the 
simple kinematic selection reported in Sec. IVI CI More- 
over, optimizing the size of the mass window, a signifi- 
cance of S/^/S + B = 1.86 is obtained in the mass win- 
dow 72 < Mjj < 110 GeV. The data needed to achieve a 
significance of 5 is reduced by a factor of two when the 
ANN selection built in this analysis is applied. 



After applying the Nout > Ncut — 0.6 cut to the 
data, the starting point of the control region (at low mjj) 
remains approximately at the same value. This essential 
feature of our ANN can be linked to the choice of having 
restricted the network training sample to events with rrijj 
values within the signal region and of not having explic- 
itly used the energy of the two leading jets in the ANN. 
In addition, the ANN cut was applied to pythia 7 + jj 
MC events to check if some discontinuity was introduced 
in the mjj spectrum between the control and the signal 
region. As expected, the mjj distribution was found to 
be very smooth over the entire mjj range. 



As far as the 



w/z 



signal distribution is concerned, af- 



ter the ANN selection, we observe no significant change 
in its Gaussian shape with the same mean and an im- 
provement of about 1% in resolution. Hence, the ANN 
has similar selection efficiency for W and Z boson events. 
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FIG. 5: Comparison of a few distributions for signal (dashed lines) and background (solid lines) events. All distributions are 
normalized to 1. The observables reported are (from left to right, top to bottom): 1) Arjjj, 2) n^ife > 3) Mj 2 /Ej 2 , 4) 5) 
f2, 6) dElf . A precise definition of these observables is provided in Table ITTT1 



E. Systematic Uncertainties on ANN Selection 
Efficiency 

Our final selection criteria are based on an ANN 
trained on simulated and real data events. Uncertain- 
ties in simulated quantities, such as jet kinematic and 
topological properties, introduce an uncertainty in the 
ANN selection efficiency. The granularity of the CDF 



detector allows an accurate determination of the direc- 
tions of jets and photons. Thus, the observables derived 
only from the directions (Arjjj, rf^ x) A$jj, maxA$j 7 , 
and min A$j 7 ) rely only upon the final state predictions 
made by the MC generator. As discussed in Sec. UH a 
good agreement on final state observables between the 
signal samples generated with pythia and madgraph 
is found; thus systematic uncertainties associated with 
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FIG. 6: Best (full square) and worst (empty square) ANN in- 
put variable combination in term of significance as a function 
of the number of input nodes for a signal efficiency of 75%. 
The fluctuations in the curves are due to small changes in 
the ANN internal parameters one has to introduce when the 
number of inputs increases. 



FIG. 8: Efficiency for the signal and the background as a 
function of the ANN output threshold Ncut- 



and ANN selection efficiency. 
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FIG. 7: The Nout normalized distributions for the signal 
(dashed line) and for the background (solid line) events. 



these variables are negligible. 



Other properties 



J 1+3 2 



Mj 2 1 E 



n ■ 



dE 3 T \ 



and 



sphericity) rely on the accuracy of the CDF detector sim- 
ulation, in particular, on the calorimeter response to par- 
ticles and track reconstruction efficiency. The calorimeter 
simulation has been extensively tuned to real data using 
isolated single tracks [l|| while track reconstruction effi- 
ciencies in data and MC are observed to be very similar. 
The dominant uncertainty on these variables comes from 
the jet energy scale. 

A change of lcr [18| in jet energy scale results in a 27% 
change in the combined 6j e t s ■ £nn signal efficiency value, 
which is assigned as total systematic uncertainty on jet 



1. Effect of Multiple pp Interactions 

In the signal sample used to train the ANN (Sec. [IT} the 
contribution of additional pp interactions (pile- up events) 
is not simulated. In this data sample the average number 
of vertices is 1.7 and more than half of the events contain 
at least one extra pp interaction. The ANN variables were 
carefully chosen to avoid any bias from pile-up events. 
The jet energies already have soft interaction contribu- 
tions subtracted (Sec. lVCj) . and only tracks coming from 
the primary vertex are considered. The only variable that 
could in principle be sensitive to additional interactions 
is the intrajet energy O since it is made up with uncor- 
rected energies. However, comparing data with single 
interaction MC 7 + jj events, a difference of less than 
4% was observed for the mean value of f2. 
To gauge the size of a possible pile-up bias in our ANN, 
we divided the data into two non overlapping sets: one 
containing events with only one reconstructed vertex, and 
the second containing events with two or more vertices. 
The ANN outputs for the two samples turned out to be 
very similar. As a further check a new ANN was built 
trained with these two samples and based on the same 
ten variables used in our analysis. With such a training 
this new ANN is built to exploit any subtle (if any) pile- 
up dependence of our input nodes and to discriminate 
events with one vertex from events with more than one. 
Similar Nout distributions (within 1%) were observed 
in the two cases, showing that pile-up events do not have 
an appreciable effect on our ANN. 
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VIII. BACKGROUND ESTIMATION AND W/Z 
PEAK SEARCH 



The search for the W/Z peak is done by subtracting 
the background contribution from the data dijet mass 
distribution. The two control regions are fitted with a 
smooth curve and interpolated inside the signal region. 
The functional form of the fit is provided by PYTHIA r )+jj 
simulated events which are best described by a simple 
exponential form f {rrijj) = e Po+Pl ' mjj . Hence, the irijj 
spectrum from the data is fitted using this function with 
Pq and Pi as free parameters. The fit is performed start- 
ing from a minimum rrijj value M m i n and excluding a 
mass window M^- < rrijj < M^ ig containing the signal 
region. For reasonable variations of these three bound- 
aries the changes in the two fit parameters were found 
to be well within their statistical uncertainties. The fit 
parameters do not show any significant change for values 
of M m i n greater than 52 GeV, while below that value 
we observe a steep increase of the fit y 2 because of the 
departure of the dijet mass shape from an exponential be- 
havior due to the trigger threshold turn-on. The fit using 
M mm = 52 GeV, M^ ig = 68 GeV, and M% g = 116 GeV 
is shown in Fig. [9] The interpolation within the signal 
region (dashed line) is our estimate of background. The 
dijet mass spectrum after the background subtraction is 
shown in Fig (TO) A consistent result for the background 
estimate was also found fitting only the high mass control 
region {rrijj > M^ g ) and extrapolating back inside the 
signal region, but at the price of a 50% larger uncertainty, 
confirming the importance of the low mass control region 
for an accurate determination of the background contri- 
bution. Since the subtracted distribution is compatible 
with zero, we are not able to identify a signal with the 
current data sample. In the next section we proceed to 
set an upper limit on the 7 + (W/Z) production with the 
W/Z boson decays into hadrons. 



IX. CROSS SECTION LIMIT CALCULATION 

To extract the signal from the data a Bayesian-based 
statistical procedure is applied. The region between 60 
and 120 GeV of the rrijj distribution is divided into 
Nbin = 15 bins, and the data events in each bin are re- 
garded as a counting experiment governed by Poisson 
statistics. The total number of events expected in the i th 
bin is Si + Bi. The number of background events Bi is 
estimated from the dijet mass distribution as described 
in Sec. IVIIII Since the stability of the control region fit 
makes the error on Bi very small, their values are held 
fixed. The number of signal events are Si = taCsi, where 
(7 is the cross section, e the total selection efficiency, C the 
integrated luminosity, and Sj the i th bin content of the 
signal dijet mass density distribution as extracted from 
the MC simulation (Fig. 0J. At first Si is held fixed as 
well. However, we show later how to take into account 
the uncertainties affecting the shape of the signal distri- 
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FIG. 9: Dijet Mass distribution of the data after the NN 
selection cut. An exponential function e Po+Pl " m " is used to 
fit the two sidebands (solid line) and the result is interpolated 
inside the signal region (dashed line). The values M m ; n =52 
GeV and [M^ g ,M p g ] = [68, 116] GeV are used to search for 
the W/Z mass peak. 
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FIG. 10: Excess of events in the data with respect to the back- 
ground prediction deduced from the sideband fits (bin errors 
do not include the background prediction uncertainties) . The 
turn-on effect can be noticed in the first two bins (they are 
not included in the background fit). No evidence of any ex- 
cess from the W/Z resonance production is found inside the 
signal region. 



bution. 

The joint probability of measuring events when 
aeCsi + Bi are expected is given by 



P{m\a,e,C) = J] 

z=l 



(aeCsi+Bi) 



In Bayesian statistics the parameters a, e and C are 
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represented by probability distributions. Before the mea- 
surement their corresponding prior density functions, 
7r(cr), 7r(e), and tt(C), summarize our a priori knowledge 
of them. Since no information on the cross section is as- 
sumed before the measurement a uniform distribution is 
chosen as its prior. In particular we define tt(ct) = if 
a < and ir(cr) = 1 if a > 0. For the efficiency and in- 
tegrated luminosity we use the estimated values eo ± Ae 
reported in Sec. I VI Bl and Cq ± AC pb -1 as reported in 
Sec. II Al Their priors are assumed to be represented by 
Gamma distributions j(x;p,, ay,) with mean /i = e 0l £ 
and width ay, = Ae, AC. The expression for the joint 
posterior probability density for (a, e, L) is provided by 
the Bayes' Theorem as: 



p(a,e,C\ni) = -^P{rii\a, e,£)7r(cr)7r(e)7r(£), 

where the normalization factor J\f constrains the inte- 
gral of p(<7, e, C\rii) to unity when integrated over all the 
parameter space. To determine the cross section we cal- 
culate the marginalized posterior probability distribution 
for cr as: 



p(a\m 



p{a, e',C'\ni)de'd£. 



However, since the jet energy scale (JES) uncertainty 
results in a change of the signal dijet mass distribution 
shape (Sec. IVII E[) . Si cannot be considered fixed and 
its dependence on JES systematics must be taken into 
account. To include this effect in the p(a\ni) computation 
a new signal density distribution is constructed moving 
the JES by one standard deviation. Its bin content is 
defined as + As.;. As a consequence the number of 
expected events is redefined as aeC(si+tAsi)+Bi, where 
the real number t parametrizes the uncertainty on the 
signal density shape. The prior density ir(t) is assumed 
to be a Gaussian distribution centered at zero and with a 
width equal to one. The posterior density for cr, including 
the new parameter t, is given by 



p(cr\m 



1 



P(n z \a, e', C, t')n(e')n(C')ir(t') de' dC dt 



As far as the cross section is concerned, this probability 
density expresses the complete summary of the measure- 
ment. Upper limits (or a central value with errors) can 
be hereby extracted from p{a\rii). The p(a\rii) distribu- 
tion was computed numerically and no local maximum 
for cr > was found. The cross section upper limit o~u m 



at 95% confidence level is computed solving the equation: 



p(a'\ni)da' = 0.95. 



It gives the value an m = 54 pb. 

X. CONCLUSIONS 

We have developed a neural network approach to iden- 
tify the dijet resonance of the W and Z boson from the 
production of events having two jets with an associated 
photon. As compared with a cut-based approach, the 
signal over background ratio improves by 163%, and the 
integrated luminosity needed for a W/Z — ► jj peak to 
emerge from the huge QCD background is reduced by a 
factor two. When applied to 184 pb -1 of data collected 
by the CDF II detector, no evidence of a W/Z — > jj 
peak is observed. The standard model prediction for 
a(pp -> Wj) X 58 (W -> qq') + <j{pp -> Zj) X 58(Z -> qq) 
is estimated to be 20.5 pb for photons with Et > 10 
GeV and \r)\ < 1.2. A 95% confidence level upper limit 
on this cross section is extracted from the data with a full 
Bayesian approach and found to be 54 pb. The technique 
employed in this analysis can be profitably extended to 
the search for small dijet resonance peaks embedded in 
large multi-jet backgrounds. 
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